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Abstract 



We discuss effects related to the fact that the final state particles of the reaction e + e~ — *> tiH 
are actually produced and they decay off mass shell. For the intermediate mass Higgs boson, 
which decays preferably into a 66-quark pair, the reaction will be observed through reactions 
with 8 fermions in the final state. Such reactions, already in the lowest order of the standard 
model, receive contributions typically from a few dozen thousands of the Feynman diagrams, 
the vast majority of which constitute background to the signal of associated production of the 
top quark pair and Higgs boson. In order to illustrate pure off mass shell effects we neglect 
the background contributions and compare the 'signal' cross section with the cross section in 
the narrow width approximation for e + e~ — > bud bpTv^ bb, which is one of possible detection 
channels of the associated production of the top quark pair and Higgs boson at a linear collider. 
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1 INTRODUCTION 



If the Higgs boson exists in Nature then it will be most probably discovered at the Large 
Hadron Collider, but the accurate study of its production and decay properties, which is crucial 
for verification of electroweak (EW) symmetry breaking mechanism, can be best performed 
in a clean environment of e + e~ collisions at the future International Linear Collider (ILC) [1] . 
The study would be of the utmost importance for establishment of the Standard Model (SM) 
or possibly some of its extensions as, e.g. the minimal supersymmetric SM (MSSM) or some of 
more general supersymmetric models as the theory of the EW interactions. The Higgs boson 
mass ma can be constrained in the framework of SM by the virtual effects it has on precision 
EW observables. The combined value of the top quark mass measured at Tevatron and the 
combined H^-boson mass [2] give a central value of m# = 85^28 GeV and an upper limit of 
166 GeV, both at 95% CL, in agreement with combined results on the direct searches for the 
Higgs boson at LEP that lead to a lower limit of 114.4 GeV at 95% CL [3J. These constraints 
indicate that the SM Higgs boson should be searched for in the mass range just above the 
lower direct search limit [I]. 

If the Higgs boson mass is below the top quark pair threshold, m# < 2m t , then the Higgs 
Yukawa coupling to the top quark g U H can be directly determined in the process of associated 
production of the top quark pair and Higgs boson [5] 

e + e- -> ttH. (1) 

The lowest order SM Feynman diagrams of reaction (pQ), with the neglect of the scalar boson 
couplings to electrons, are shown in Fig. [TJ As the contribution of the Higgs boson emission off 




Figure 1: Feynman diagrams of reaction (CQ) to the lowest order of SM with the neglect of the 
scalar boson couplings to electrons. 

the virtual Z-boson line, which is represented by the diagram in Fig. [lb, is small with respect 
to the Higgsstrahlung off the top quark line illustrated in Fig. [T^l and [lb, the SM lowest order 
cross section of reaction ([I]) becomes practically proportional to g^ tH . This fact makes reaction 
([T]) so sensitive to the Higgs Yukawa coupling to the top quark. Because of its numerical value 
close to 1, precise determination of g tt u m &y play an indispensable role in our understanding 
of the whole mass generation mechanism of SM. 

Because of their large decay widths, the t- and t-quark of reaction (CQ) almost immediately 
decay into bW + and bW~, respectively, the H^-bosons subsequently decay into 2 fermions 
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each and the Higgs boson, if it is lighter than 140 GeV, decays dominantly into a 66-quark 
pair. Thus reaction ([I]) will be actually detected at the ILC as a reaction of the form 

e+e- - bfj[bf 2 ftbb, (2) 

where fi, = v e , ^> v T , u, c and f[, fi = e~, r~, d, s. The three possible detection channels 
of (j2J), which correspond to the decay modes of the TV-bosons resulting from decays of the t- 
and t-quark in the 'signal' diagrams of reaction (jTJ), are 

1. hadronic channel: eight jets (38 %), 

2. semileptonic channel: lepton and six jets (37 %), 

3. leptonic channel: two leptons and four jets (25 %). 

In addition to the analysis in ref. [5], reaction ([ID has been extensively studied in literature. 
The QCD radiative corrections to (JTJ) were calculated in [6], 0(a) EW corrections were calcu- 
lated in [7] and full 0(a) EW and 0(a s ) QCD corrections were studied in [S]. In [5], process 
(JTJ) was considered in the kinematic region where the Higgs boson energy is close to its max- 
imal energy and hence the next-to-leading-logarithmic corrections to the Higgs boson energy 
distribution can be computed within the nonrelativistic effective theory. Processes of the form 
e + e~ — > bbbbW + W~ — > bbbbl^uiqq' accounting for the signal of associated Higgs boson and 
top quark pair production, as well as several irreducible background reactions, were studied in 
[TU] and EW contributions to the leptonic and semileptonic reactions (j2j) have been computed 
in [11]. Moreover, feasibility of the measurement of the Higgs-top Yukawa coupling at the ILC 
in reaction flTJ was discussed in [12]. 

Reactions (J2]), with 8 fermions in the final state, receive contributions typically from a few 
dozen thousands of the Feynman diagrams already in the lowest order of SM. Most of the 
diagrams constitute the 'non signal' background to the reaction of the associated on shell top 
quark pair and Higgs boson production and their subsequent decay 

e + e - uH - bhf[bhf 2 b\ (3) 

with the same final state as that of (l2l). 

In the present note, we will look into pure off shell effects that are related to the fact that t, i 
and H of ([3]) are actually produced and they decay as off mass shell particles. To illustrate these 
effects, we will calculate the cross section of one selected reaction ([2]) while keeping only the 
'signal' Feynman diagrams and compare it with the cross in the narrow width approximation 
(NWA) for the top, antitop and Higgs in the energy range that could be relevant for the ILC. 
We will also check whether or not the off shell effects change the prediction of [5] that the 
cross section of ([T]) is dominated by the Higgs boson emission off t and t. We realize that 
neglecting an overwhelming number of the 'non signal' Feynman diagrams may be too crude 
an approximation, in particular in the energy range not far above the ttH threshold. However, 
our simplified approach has the advantage that it fully takes into account spin correlations 
that are of great importance in the context of top quark physics [13] . Taking into account the 
spin correlations would in particular increase the sensitivity to new physics effects, such as an 
anomalous Wtb couplings [H) . especially if the beam polarization is provided at the ILC. 
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2 OUTLINE OF THE CALCULATION 



We will concentrate on one selected semileptonic channel of reaction (j2J), namely 

e + Oi) e~(p 2 ) -> b(p 3 ) u(p 4 ) d(p 5 )b(p 6 ) fi'(p 7 ) u^ps) b(p 9 )b(p 10 ), (4) 

where the particle four momenta have been indicated in parentheses. In the unitary gauge, 
reaction (j3J) receives contributions from 56550 Feynman diagrams already in the lowest order 
of SM. The corresponding amplitudes can be generated with a Fortran 90 program carlomat 
written by one of the present authors [15] . If we neglect the Higgs boson couplings to fermions 
lighter than a 6-quark then the number of diagrams is reduced to 26816. However, as already 
stated in Section 1, in order to illustrate the size of the pure off mass shell effect we will restrict 
ourselves to the lowest order 'signal' Feynman diagrams shown in Fig. [2j The corresponding 




Figure 2: 'Signal' Feynman diagrams of reaction (J4]). 



amplitudes are calculated with the helicity amplitude method that had been described in 
detail in [16] and computed with the use of program libraries described in [TT] . 

With that restricted numbers of amplitudes, the actual challenge of the computation of the 
total cross section of (jl]) is a right choice of parametrization of the 20-fold phase space integral 
which should lead to reliable results without using vast numbers of calls to the integrand. The 
3 different phase space parametrizations which we use are the following 

d 2 °LipS = (27r) 20 dPi?2 ( s j s 345910> s 67s) dPS 2 (S345910, $345) s 91o) 

x dPS 2 (s 3 45, ml, s 45 ) dPS 2 [s 678 , ml, s 7 fj 

x dPS 2 (s 4 5, ml, mfj dPS 2 (s 78 , m 7 , mfj dPS 2 (s 9 i , ml, m\^j 

X dS34591ods 6 7 8 dS345dS45ds 78 dS9io, (5) 

d 20 Lips = (27vy 20 dPS 2 (s, s U5 , s 6789 io) dPS 2 (s 345 , m 2 3 , s 45 ) 
x dPS 2 (s 45 , ml, mfj dPS 2 (s 6789 i , s 678 , s 910 ) 
x dPS 2 (s 678 , ml, s 78 ) dPS 2 (s 78 , m 7 , mfj dPS 2 (sg 10 , mj, mffj 

X dS345dS67891odS45dS678ds 78 dSgio, (6) 

d 20 Lips = (2ny 20 dPS 2 (s,s M5e78 ,s 9W )dPS 2 (s U5678 ,s U5 ,s 678 ) 

dPS 2 (s 34 5, ml, s 45 ) dPS 2 (s 45 , m\, mfj dPS 2 (s 678 , m\, s 78 j 



x 
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x dPS 2 (s 78 , to 2 ., mgj dPS 2 [s 910 , mj, m\ Q 

x ds345678dS345dS678dS45dS7gdS9io. (7) 

In Eqs. s = {pi+p 2 ) 2 , = (Pi + Pj + -) 2 , h J = 3, ...,10, and dPS 2 (g 2 , ft 2 , q 2 ) is a 

two particle (subsystem) phase space element defined by 

dPS 2 (q\ ft 2 , ft 2 ) = (q - ft - ft) |||| = ^d^, (8) 

where ft is the momentum and tt% is the solid angle of one of the particles (subsystems) in 
the relative centre of mass system, ft + q 2 = 0. 

Parametrizations (j5J), (EI) and (JTj) have been chosen in such a way that invariants Sy... corre- 
spond to virtualities of propagators of the gauge bosons, Higgs boson and/or top quarks in the 
diagrams of Fig. [2j Possible poles in the propagators of unstable particles are regularized with 
the constant particle widths T a which are introduced through the complex mass parameters 
M 2 by making the substitution 

m 2 a -> M 2 = m 2 a -im a T a , a = Z,W,H,t. (9) 

In order to reduce variance of the MC integration invariants s^... that are related to the 
resonating propagators of W, H and t are obtained from the random variables uniformly 
distributed in the interval [0, 1] by performing mappings smoothing out their Breit-Wigner 
distributions. Denote the lower and upper physical limit of s^-... by s^ m and s™ ax , respectively, 
and the uniform random variable by x G [0, 1], then the mapping is given by 

Si,-... = T a m a tan ( ^ a ™ a x + Xq) + m 2 a , (10) 



N 



with 



-.mm 



N = /, m ax m2 V and x o = ~ arctan — J — . (11) 

arctan ( _ Xq ) \ T a m a / 

Mapping (TX0|) is not performed for S345910, S678910 and S345678 hi the propagators of t in Fig. [2^, 
t in Fig. [2b and Z in Fig. [2k, respectively, which are far away from resonance in the centre 
of mass system energy (CMS) range considered. Those three invariants, as well as all the 
remaining integration variables of Eqs. (JSD — (JHJ) , which we symbolically denote by yj G [aj,bj], 
are obtained from the uniform random variables Xj G [0, 1] with a simple linear transformation 

y j = (b j -a j )x j + a j . (12) 



Weights Wi, i = 1,2,3, with which each of the 3 phase space parametrizations (jSI)-© con- 
tributes to the total cross section are determined in the initial scanning run performed with 
all the initial weights equal to 1/3. They are calculated as the ratios 

3 

w i = a i /^a j , i = 1, 2, 3, (13) 



5 



where <7j, i — 1,2, 3, denotes the cross section obtained in the initial scan with phase parametriza- 
tion (j3J), (jBJ), (JTj), respectively. The final result for the total cross section a of (j4j) in the 'signal' 
approximation is then calculated as the weighted average 

3 

a = J2 w i (T j ( 14 ) 

3=1 

with aj, j = 1, 2, 3 being the cross section computed with phase parametrization ((Sj), (J6]), (JTJ), 
respectively. 

The cross section of (j3j) in the NWA, is defined in the following way 

o-nwa = o"(e + e -> ttiij x — x £ x — , (15) 

1 w >-W J- H 

where a(e + e~ — ► tDf ) denotes the total cross section of ([1]) and we have assumed 

Tt-tbw+ Ff-fiw- 



1. 



The 5-dimensional numerical integration that is necessary in order to compute a(e + e~ — > tiH) 
does not require the multi-channel MC approach and can be performed with the single phase 
space parametrization given by 

d 5 Lips = (2n)- 5 dPS 2 (s,s t t,m 2 H )dPS 2 (su,m 2 t ,m*)dsu, (16) 



where s t i = {pt +Pt) 2 , with p t and pi being the four momenta of the on shell t and i in (JTJ 
and the two particle phase space element dPS 2 defined in (jSj). 

The NWA formula of Eq. ( Tl5|) is obtained from the signal cross section of reaction (j4]) by 
making the following substitution for the resonance factors D a (g 2 ), a = t, W, corresponding 
to each of the resonating propagators of the top quark and IV-boson in the Feynman diagrams 
of Fig. m 



Da[q 



<? - m a) + (ra a r a ) 



K a 6(q 2 -m 2 a ), (17) 



with normalization factors K a = J dq 2 D a (q 2 ) — tt/ (m a T a ). It is just the Dirac delta function 

— oo 

on the right hand side of (1171) which causes the signal cross section of (jl]) to take the factorized 



form of (jl5p . Obviously, substitution (ITT)) makes sense only if the width of unstable particle 
r a is much smaller than its mass m a . The relative error associated to it is usually estimated 
as being of O (T a /m a ), as required by dimensional analysis. With the values of the top quark 
and H^-boson masses and widths given in Eqs. (IT51) . (12111) and (|2"Tj) we obtain uncertainties 
of 0.9% and 2.5% for each of the resonating top quark and jy-boson propagators. Hence we 
may expect a discrepancy between the signal and NWA cross sections of the order of a few 
per cent. This expectations will be confirmed by the actual numerical results which are given 
in the next section. For a more detailed discussion of limitations of the NWA the reader is 
referred to 
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3 NUMERICAL RESULTS 



The numerical results presented in this section have been obtained with the following set of 
initial physical parameters. We have chosen the Fermi coupling and fine structure constant in 
the Thomson limit 

G M = 1.16639 x 1CT 5 GeV -2 , a = 1/137.03599911, (18) 

as well as the W- and Z-boson masses 

m w = 80.419 GeV, m z = 91.1882 GeV, (19) 

as the EW input parameters. The top quark mass and the external fermion masses of reaction 
(j3J) are the following: 

mt = 174.3 GeV, m b = 4.8 GeV, m u = 5 MeV, m d = 10 MeV, (20) 
m e = 0.51099892 MeV, m M = 105.6583 MeV. 

The value Higgs boson mass is assumed at m# = 130 GeV. The widths of unstable particles 
that are introduced through substitution (Q are calculated to the lowest order of SM. This 
results in 

T t = 1.531 GeV, T w = 2.048 GeV, T H = 8.065 MeV, (21) 

for the top quark, W-boson and Higgs boson widths, respectively. The actual value of the 
Z-boson width is not relevant for our calculation, as the Z-boson propagator is far off its mass 
shell in all the diagrams of Fig. [2J 

The total cross section of reaction (j3J) calculated with the 'signal' diagrams of Fig. [2] is com- 
pared with the cross section in the NWA in Table [H We see that the off mass shell effect 
is of the order of 3% for y/s = 500 — 800 GeV, i.e. close to the threshold of the associated 
top quark pair and Higgs boson production, then it decreases and becomes negative, reaching 
an absolute value of about 5% at \fs = 2 TeV. This is illustrated in Fig. [3j where both the 
'signal' and NWA cross sections are plotted as functions of the CMS energy. 

The off shell effect is caused overwhelmingly by the nonzero widths of the top quarks and 
VT-bosons, as for the assumed value of the Higgs boson mass, m# = 130 GeV, the Higgs 
boson is very narrow, with a decay width of a few MeV. Actually, the very small lowest order 
value of 8 MeV for T H in Eq. (I2T]) will be further substantially reduced if the QCD radiative 
corrections are taken into account, see, e.g. [19] and |20j . 

It is interesting to see to which extent the off shell effects may change the prediction of [5j that 
the cross section of (CQ) is dominated by the Higgs boson emission off the t- and t-quark. To 
this end, in Table [21 we present the lowest order 'signal' cross section of PJ o" signa i and the cross 
section crf° g ^ z that has been calculated with the diagrams of Fig. [2^ and [2b, i.e. without the 
Higgsstrahlung off the Z-boson line represented by the diagram of Fig. [2b. The corresponding 
results for the cross section in the NWA are given in Table [2] as onwa and o"nwa ZZ - Let us 
compare the relative differences shown in Table [2] as 5\ and 62 defined by 

$1 = (^SgrSf 2 — ^signal) /^signal, #2 = (cnWA^ — 0"NWa) /CNWA- (22) 
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[GeV] 


^signal [ab] 


ffNWA [ab] 


5[%] 


500 


3.805(11) 


3.923(1) 


3.1 


800 


58.33(7) 


60.07(3) 


3.0 


1000 


51.79(7) 


52.56(3) 


1.5 


1200 


42.99(6) 


42.96(3) 


-0.1 


2000 


21.90(11) 


20.76(2) 


-5.2 



Table 1: Total cross sections of reaction (TjJ: the 'signal' cross section o" signa i, the cross section 
in NWA ctnwa and their relative difference 6 = (onwa — ^signal) /^signal- The numbers in 
parenthesis show the uncertainty of the last decimals. 

They quantify to which extent the Feynman diagram of Fig. [2fc spoils proportionality of the 
cross section of (Jl]) to g^ H , which makes the measurement of the top-antitop-Higgs Yukawa 
coupling more difficult. While 8\ takes into account the fact that t, i and H are produced and 
decay off mass shell, 62 has been calculated assuming that they are on mass shell particles. 
The lowest order 'signal' cross section of (jl]) o" S i gna i and the 'signal' cross section without the 
diagram of Fig. [2fc cr^g^f 2 are plotted as functions of the CMS energy on the left hand side 
of Fig. HI The plots on the right hand side of Fig. H] show the corresponding cross sections in 
the NWA. Both from Table [2] and Fig. HJ we see that the Higgsstrahlung off the Z-boson line 
spoils the proportionality of the total cross section of the associated tiH production to g\ m 
almost in the same way, independently of whether the off shell effects are taken into account 
or not. 



4 Summary and outlook 

We have looked at a role that the off shell effects may play in reaction ([T]) of the associated 
production of the top quark pair and Higgs boson at the ILC. We have illustrated these effects 
for a semileptonic reaction (T4]), which is one of the detection channels of at the ILC, by 
comparing the 'signal' cross section that has been calculated by performing 20-fold integration 
of the squared matrix element while keeping only the 'signal' Feynman diagrams with the cross 
section in the NWA. The off shell effects are typically of the order of a few per cent for the 
CMS energies in the range from 500 GeV to 2 TeV, in accordance with the expectation based 
on the discussion of the quality of the NWA in the end of Section 2. We have also shown 
that the off shell effects do not affect much the prediction of [5j that the cross section of ([T|) 
is dominated by the Higgs boson emission off t and i, which makes it an attractive tool for 
determination of the Higgs-top Yukawa coupling. 

The presented approach is very simplified, as a lot of the 'non signal' background Feynman 
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Figure 3: Total cross sections of reaction as functions of the CMS energy: the 'signal' cross 
section cr signa i (solid line) and the cross section in NWA cxnwa (dashed line). 




•\A (GeV) y-s (GeV) 

Figure 4: The lowest order cross sections of (j3J) as functions of the CMS energy. On the 
left hand side: the 'signal' cross section (solid line) and the 'signal' cross section without the 
diagram of Fig. (dashed line); on the right hand side: the corresponding cross sections in 
the NWA. 
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1 — 

[GeV] 


""signal 

[ab] 


nn HZZ 
signal 

[ab] 


Si 
[%) 


CNWA 

[ab] 


nn HZZ 

°NWA 

[ab] 


c 

$2 
[%) 


500 


3.805(11) 


3.775(10) 


-0.8 


3.923(1) 


3.886(1) 


-0.9 


800 


58.33(7) 


55.84(6) 


-4.3 


60.07(3) 


57.51(3) 


-4.3 


1000 


51.79(7) 


48.52(6) 


-6.3 


52.56(3) 


49.23(3) 


-6.3 


1200 


42.99(6) 


39.50(7) 


-8.1 


42.96(3) 


39.41(3) 


-8.3 


2000 


21.90(11) 


19.18(12) 


-12.4 


20.76(2) 


18.03(2) 


-13.2 



Table 2: Total cross sections of reaction (TJJ: the 'signal' cross section o" signa i, the 'signal' 
cross section calculated without the diagram of Fig. |2fc o-f° g ^ u , the cross section in the NWA 
c"nwa, the cross section in the NWA calculated without the diagram of Fig. [2b ^nwa 11 an< ^ ^ e 
relative differences Si and S2 of ( 1221) . The numbers in parenthesis show the uncertainty of the 
last decimals. 

diagrams have been neglected in the calculation of the lowest order cross section of (j4j), but it 
has the advantage that it fully takes into account spin correlations that are of great importance 
in the context of top quark physics. Further work is needed in order to take into account a 
complete set of the lowest Feynman diagrams and include leading radiative corrections for 
reactions (T5]). 
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